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Abstract 



We study a flavour- violating four-fermion interaction in the Lifshitz context, in 3+1 di- 
mensions and with a critical exponent z = 3. This model is renormalizable, and features 
dynamical mass generation, as well as asymptotic freedom. At one-loop, it is only loga- 
rithmically divergent, but the superficial degree of divergence of the two-point functions 
is 3. We calculate the two-loop corrections to the propagators, and show that, at this 
^ ■ order, the Lorentz-violating corrections to the IR dispersion relation are quadratic in the 
cut off. Furthermore, these corrections are too important to represent a physical effect. As 
a consequence, the predictive power of the model in terms of Lorentz-violating effects in 
! the propagation of particles is limited. 

O 



1 Introduction 

Lifshitz-type theories, where time and space have different mass dimensions and therefore 
violate Lorentz invariance, have attracted attention in recent years, motivated essentially 
by the possibility of defining new renormalizable interactions. This is because of an im- 
provement in the convergence of loop integrals, due to higher order space derivatives, which 
is achieved without the introduction of ghost degrees of freedom, since the order of time 
derivatives remains minimal. A review of Lifshitz theories can be found in [1] for quantum 
field theories in particle physics and in [2] for the Horava-Lifshitz alternative approach 
to gravity. An example of a new renormalizable interaction in the Lifshitz context is the 
Liouville interaction in 3+1 dimensions [3], where the exponential potential for the scalar 
field is a relevant interaction because the field is dimensionless if the critical exponent is 
z = 3. 
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We consider a Lifshitz-type four-fermion interaction model, which, in d = 3 space di- 
mensions and for an anisotropic scaling z = 3, is renormalizable [4]. Such models have 
been studied in [5], where two fundamental properties were shown: dynamical mass genera- 
tion and asymptotic freedom. One motivation for these theories, besides renormalizability, 
is the apparent improvement of quantum corrections, since the one-loop graphs are only 
logarithmically divergent, instead of quadratically in the Lorentz case. But the overall 
superficial degree of divergence of the graphs of the theory is actually ou = 6 — 3E/2, where 
E is the number of external lines. If one considers the propagator (E = 2), the corre- 
sponding corrections have a superficial degree of divergence equal to 3. The coefficient of 
the cubic divergence may cancel for some graphs, but we calculate here a two-loop graph 
which shows that the divergence in the model is at least quadratic. Therefore, although 
renormalizable, this Lifshitz model still contains "large" divergences. 

Our model features two massless fermion flavours, coupled with four-fermion interac- 
tions which do not respect flavour symmetry. After showing the occurrence of dynamical 
flavour oscillations in this model, we calculate the modified dispersion relations for these 
two fermions, arising from quantum fluctuations. Classically, all fermions have the same 
dispersion relations, with higher order powers of the space momentum p, rescaled by a large 
mass M, which represents the crossover scale between the Lifshitz and Lorentz regimes. 
These dispersion relations coincide with the expected Lorentz-invariant one in the infrared 
(IR) regime \p\ « M. Taking into account quantum corrections, though, modifies this IR 
limit: it is known in Lifshitz-type studies that different species of particles see different ef- 
fective light cones [6]. Since our model breaks flavour symmetry, the dressed IR dispersion 
relations are different from the Lorentz-invariant one, and we show that the corresponding 
corrections are quadratically divergent. Furthermore, as we will see, these corrections are 
too important to represent any physical effect. As a consequence, a proper treatment of 
the model would consist in defining counterterms to absorb these divergences, and no pre- 
diction can be made as far as Lorentz-violating propagation is concerned (a logarithmic 
divergence could lead to a "realistic" energy dependent effective maximmum speed). 

The next section derives the dynamical masses for the model, including the mass mix- 
ing terms necessary for flavour oscillations. From our study of dynamically induced flavour 
oscillations, and taking into account experimental data on neutrino oscillations, we de- 
rive values for the coupling constants of our model, which, as expected, are perturbative. 
Although neutrinos are not Dirac fermions, the corresponding experimental constraints 
give a good order of magnitude for the parameters in our model. Section 3 shows the 
asymptotic freedom of the interaction, based on a one-loop calculation. The four-point 
function has a vanishing superficial degree of divergence, such that higher order correc- 
tions cannot change the sign of these beta functions. The effective IR dispersion relations, 
dressed by quantum fluctuations, are derived in section 4. For this, we need to go to two 
loops, since the one- loop correction to the fermion propagators is momentum-independent. 
Detailed calculations are given in the Appendix, where we perform part of the integration 
analytically and then integrate the rest numerically. 
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2 Flavour oscillations 



We describe here, in the Lifshitz context, how flavour oscillations can arise dynamically 
from flavour-mixing interactions between two massless bare fermions, as suggested in [7]. 
From our expressions for the dynamical masses, together with experimental data, we find 
phenomenologically realistic values for the coupling constants of our model. 
Oscillations of massless neutrinos are studied in [8], where neutrinos are considered open 
systems, interacting with an environment. Such oscillations have also been studied in [9], 
in the framework of Lorentz- violating models, involving non-vanishing vacuum expectation 
values for vectors and tensors. Whilst these studies have been questioned by phenomeno- 
logical constraints [10], our present model, based on anisotropic space time and higher 
order space derivatives, is not excluded. 

Flavour oscillations were also related to superluminality in [11], where it is shown that, if 
superluminality is due to a tachyonic mode, the latter can be stabilised by flavour mixing. 
Finally, in [12], superluminal effects are related to the extension of a single neutrino wave 
function, where the oscillation mechanism plays a role in the uncertainty of the neutrino 
position. 

2.1 Flavour symmetry violating 4-fermion interactions 

We work in the z = 3 Lifshitz context, in d = 3 space dimensions. We consider two flavours 
of massless Dirac fermions ipi,ip2, and the free action 

Sfree = J dtdx^fj^ - ^ a (M 2 - A)(id ■ 7)^ a = 1, 2 , (1) 

where [M\ = 1 and [ip a ] =3/2, and a dot over a field represents a time derivative. For 
the dispersion relations to be consistent in the IR (see eq.(4) below), one can consider M 
typically of the order of a Grand Unified Theory scale (GUT), although we will show that 
our results only slightly depend on the actual value of M. We introduce the following 
renormalizable, flavour-violating and attractive 4-fermion interactions 

S in t = j ' dtdx(g 1 i(> 1 ip 1 + 02^2^2 + M^Mta + » ( 2 ) 

where the coupling constants g±,g2,h are dimensionless. As shown in [5], this kind of 
model exhibits dynamical mass generation, which can be seen only with a non-perturbative 
approach, as will be shown in the next section. Taking into account the dynamical masses, 
but ignoring quantum corrections to the kinetic terms, the dispersion relations are of the 
form 

^ 2 = <„ + (M 2 + P 2 )V, a = 1,2, (3) 
which, after the rescaling u = M 2 u, leads to 



where rhdyn = m d yn /M 2 . One can see then that Lorentz-like kinematics are recovered in 
the IR regime p 2 « M 2 , as expected in the framework of Lifshitz models. After the 
rescaling t = t/M 2 



<< M 2 , as expected in the framework of Lifshitz models. 
:1 the action reads 



= J dtdx ^ J $ip a + r„-^7( /V • 7)< '„ 



+ 



+ ||^2^2 + ^Ol^2 + Ml] 



(5) 



where we can see that the four fermion couplings (g a /M) 2 , g a h/M 2 and (h/M) 2 are very 
small compared to the Fermi coupling ~ 1CT 5 GeV~ 2 , if M is of the order of a GUT 
scale, or even several orders of magnitude smaller, and g a , h are perturbative. Finally, note 
that, for Large Hadron Collider energies up to few TeVs, the classical Lifshitz corrections 
p 4 /M 2 and p 6 /M 4 in the dispersion relation (4) are not detectable, if M is of the order of 
a GUT scale. For this reason, if one wishes to describe measurable non-relativistic effects 
in the Lifshitz context, these should be sought in quantum corrections to the IR dispersion 
relation. 



2.2 Superficial degree of divergence 

It is interesting to note that, although this Lifshitz model has only logarithmic divergences 
at one-loop, quantum corrections actually do not "behave better" than those in the Lorentz- 
invariant 4 theory, since the superficial degree of divergence of the propagator is 3. 
To show this, we calculate via the usual approach the degree of divergence a; of a graph 
with E external lines. Each loop gives an integration measure dp d 3 p, which has mass 
dimension 6, and each propagator has mass dimension -3. For a graph with / internal 
lines and L loops, the superficial degree of divergence is therefore u = 6L — 31. As usual, 
because of momentum conservation, we also have L — I — n + 1, where n is the number 
of vertices of the graph. Finally, since we have 4-leg vertices, we also have the relation 
4n = E + 21. Taking into account these constraints, we find u = 6 — 3E/2. 
From this result, we see that the four-point function is at most logarithmically divergent, 
but the propagator has a superficial degree of divergence equal to 3, although the one-loop 
mass corrections are logarithmically divergent only, as we show in the next subsection. 



2.3 Dynamical generation of masses 

We now calculate the dynamical masses generated by the interaction (2). For this, we 
introduce the auxiliary scalar field (f> to express the interaction as 

exp (iS int ) = J V[4>] exp(iS' ) , 
with = / dtdx (-0 2 + 20(5-1^1^1 + 92^2^1 + ^(^1^2 + ^2^1)) , (6) 
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and then calculate the effective potential for 0=constant as 

exp (iVV r e// (0)) = J V[il) 1 ,tl) 1 , ^2,^2} e*P(iS free + iS<p) , (7) 

where V is the space time volume. This integration can be done exactly, since Sf ree + is 
quadratic in fermion fields, and leads to an effective potential for <fi. From the dispersion 
relation (4), one can see that a non-trivial minimum <p m i n for this effective potential will 
give the flavour mixing mass matrix 

(# l)< < 8 > 

which leads to the rescaled masses 

(fill P> \ _ n ^min ( 9l h \ , v 

{ fi m 2 ) A { h g 2 ) ■ W 
As a consequence, the mass eigenstates are 

m± = (^i + S2 ± (51 - ^)v/l + tan 2 (2^)) , (10) 

where the mixing angle is defined by 

tan(20) = — ^— • (11) 
9i ~ 92 

With the auxiliary field, the Lagrangian can then be written in the form tyO^, where 



»= £ . (») 



and the operator is 



/ i^d -(M 2 -A)(id^) + 2g 1 cf ) 2h4 \ ( , 

V v 2h<f) * 7 °<9o-(M 2 - A)^-7) + 2^0 / ' 1 J 

Integration over the fermions then gives the following effective potential for <fi (where the 
Euclidean metric is used for the loop momentum) 



1 r duo dp 



2 J 27r(27r) 3 



VefAt) = ^--J tFtt^^ ^ + (M 2 +p 2 )V (14) 
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+ 40 2 c 2 + (M 2 + p 2 ) V fo? + ^ 2 + 2/> 2 ) + 16(^ - ** 2 ) V 



A derivative with respect to gives 



d(f) r W 27r(27r) 3 (w 2 + C + )(a; 2 + C_) ' 



where 

A = A{gl + gl + 2h 2 ) (16) 
B = 4(M 2 +p 2 ) 2 p 2 (g 2 1+ g 2 2 + 2h 2 ) + 320 2 ( 9l g 2 -h 2 ) 2 

C± = (M 2 + p 2 ) V + 20 2 [(^ 2 + ^ 2 + 2/> 2 ) ± y^ 2 - glf + + g 2 f . 

(17) 

The integration over frequencies u leads to 

dV eff 1 /" A 2 , ( B A 



do =2 *- +wk " dp [cWcCTcWc; + TcTWcl) (18) 

where A is the UV cut off, assumed to be large compared to M. A non-trivial minimum 
<t>min 7^ for this effective potential is solution of the equation 



The dominant contribution of these logarithmically divergent integrals comes form p — > A, 
so we can therefore approximate 

C_ ~ C+ ~ p 6 + ^0 2 and B~Ap 6 , (20) 



such that 



The gap equation (19) then gives 



A 3 / -6tt 2 \ , 

<Pmin ^ / = — exp -2— — , (22) 

^g\ + ^ 2 2 + 2h 2 \g 2 1 +g 2 2 + 2h 2 J 



and the rescaled masses (9) are 



2g a A 3 / -67T 2 \ 

m - - r=„ = = T7^ ex P L 2 , ,2 , ou ( 23 ) 



+ ^ + 2^ ^ 2 W + sl + 2/* 2 , 



2h A^ / -6tt 2 \ 

" " y/g} + gl + 2h* M 2 ^ [g 2 + , 2 2 + 2tf J " (24) 

As expected, these masses are not analytical in the coupling constants and could not 
have been obtained with a perturbative expansion. Similar results have been obtained in 
the context of magnetic catalysis [14], based on the Schwinger-Dyson approach, and also 
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for Lorentz- violating extensions of QED [15]. Neither of these studies however feature 
anisotropic space time studied herein. 

Finally, we note that the approach adopted here, based on the effective potential for the 
auxiliary field 0, is in principle valid for a large number of flavours. Indeed, this auxil- 
iary field depends on space and time, and its fluctuations around the minimum 4> m in induce 
new fermion interactions. For N fermion flavours though, these fluctuations are suppressed 
by 1/N, which justifies the approach. In our case, N = 2 is not "large", but the corre- 
sponding order of magnitude for the dynamical masses is sufficient for a suitably accurate 
determination of the coupling constants <?i,<?2, as explained in the next subsection. 



2.4 Experimental constraints 

From the expressions (10), we obtain the following difference of mass eigenstates squared 

4 



9 2 i - gl 



cos(2#) g\ + gl + tan 2 (2^)(( ?1 - g 2 f/2 

A 6 / -127T 2 



x — — exp 



M 4 r\gl + gl + t^(29)( gi -g 2 y/2, 



(25) 



Experimental constraints are [16] 



Am 2 12 



7.59(7.22 - 8.03) x 10~ 5 (eV) : 



sin 2 0i2 = 0.318(0.29 - 0.36) , 



(26) 



and we plot in fig.(l), from the expression (25), the set of points in the plane g±, g 2 which are 
allowed, given the experimental constraints (26). We consider A ~ 10 19 GeV, corresponding 
to the Plank mass. An important property is that the result is hardly sensitive to the value 
of the mass scales M: because of the exponential dependence in eq.(25), an increase of 
several orders of magnitude in M leads to an increase of few percent only for the couplings 
g a , as shown in the following table. Considering the situation where h « 1, such that 
g 1 ~ g 2 , according to eq.(ll), the approximate common value for the coupling constants 
as a function of the ratio M/A is then: 



M/A 


1Q -16 


lO" 15 


10 -14 


lO" 13 


IO" 12 


io- 11 


9i - 92 


0.46 


0.47 


0.48 


0.48 


0.49 


0.50 


M/A 


io- lu 


io- y 


10~ 8 


io- 7 


IO" 6 


10~ 5 


9i - 92 


0.51 


0.53 


0.54 


0.55 


0.56 


0.58 



Table 1: Coupling constants for different mass scales M, when h « 1 and A = 10 19 GeV. 
On fig.(l), the thin line represents the set of points satisfying the constraint 



In 



^^experimental \ 
, '""calculated 



) 



< 1 , 



(27) 



-0.5 0.0 0.5 

gi 

Figure 1: Values of g\ (x-axis) and g 2 (y-axis) allowed by experimental constraints, for 
Af/A = 10 -11 . Negative values are allowed, since the physical quantities depend on the 
square of the coupling constants only. Points where g 1 = g 2 are strictly speaking not 
allowed, since at these points Am 2 = 0. However, the resulting logarithmic singularity is 
very localised in the parameter space, such that we can safely chose g\ and g 2 perturbatively 
close to each other. 

and the thick line represents the set of points such that the largest mass eigenvalue is 
between 10~ 3 and 1 eV. We see that the coupling constants appearing in the theory are 
then of the order g\ ~ 0.25, and can be considered perturbative. 

3 Asymptotic freedom 

We now calculate the one- loop coupling constants, for h « 1, and we show that the 
theory is asymptotically free. For simplicity, we set h = but still keep g\ ^ g 2 . 
The bare interaction can be expressed as 

9 2 i (^i) 2 + G V^i^2 + gl (^ 2 ) 2 , G = 2 9l g 2 , (28) 

and the dressed interaction is of the form 

{gl + ^i 2 )(?A) 2 + (G + 5G)^^ 2 + {g 2 2 + 5gl)^ 2 ) 2 . (29) 

Note that no symmetry imposes any relation between SG and 5gf,5g 2 '- the interaction 
ipi4'i' l p2' l p2 is dressed independently of the interactions (ip^i) 2 and (ip 2 ip2) 2 - 
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3.1 One-loop Fermi coupling 

If one denotes 



D a (co,p) = u; 2 -(M 2 +p 2 )Y-ml, (30) 



the generic graph for the one-loop corrections is 

dudp iN a (u,p)iN b (u,p) 
(2^p D a (u,p)D b (u,p) ' 



Lab 



(31) 



When A >> M, we obtain 

r dudp ( ml/ (ml - m|) m 3 b /(m 3 a - m|) \ 

ab 7 (2tt) 4 I a; 2 - (M 2 + V 2 W - m* uj 2 - (M 2 + v 2 ) 2 P 2 -ml) 1 ' 




47r 2 (m^ — ml) 

i 

127r 2 (m^ — mf) Jo 
i 

47r 2 (m^ — m\ 

The integral (32) diverges logarithmically, unlike the Lorentz symmetric case where it 
diverges quadratically. Note that, when m b — > m a , the previous result is regular and leads 
to 

^-T^lnf-V (33) 



4vr 2 

In order to calculate the number of graphs (32) contributing to the coupling corrections, 
we introduce the auxiliary field a and write the four-fermion interactions in the form 

- l -o 2 + av^iV^i + 52^2^2) • (34) 

The scalar a does not propagate, but is described by a fictitious propagator, which carries 
a factor i. This propagator has to be understood in the limit where it shrinks to a point, 
leading to the fermion loops given by the expressions (32). The two vertices corresponding 
to the effective Yukawa interactions are i\/2gi and i\f2g 2 . 

The graphs corresponding to the four-point function are represented in fig. (2), in terms of 
the equivalent Yukawa interaction (34), where the last two graphs do not contribute to the 
four-fermion beta functions. Indeed, the general structure of the four point function is 







^V^V 0) = Al <g) 1 + B t l <g) a 1 + djcx 1 ® a 3 , (35) 
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Figure 2: One-loop graphs involving the auxiliary scalar field, which contribute to the four- 
point function. Solid lines represent fermions and dashed lines represent the scalar. Only 
the first two diagrams, where the fermion lines cross an odd number of vertices, contribute 
to the four-fermion beta functions. The first graph corresponds to the insertion of a scalar 
self-energy, for each flavour, and involves a factor -4 for the trace over Dirac indices. The 
second graph has two contributions: one for each insertion of a vertex correction. 

where the Dirac indices are omitted and the tensorial product allows for the two in and 
two out states. A is the only quantity contributing to the coupling constant, since the 
corresponding term has no Dirac structure. The four-point function contains one divergence 
only, which is logarithmic, such that the divergent graphs are obtained only from the highest 
power of momentum in the numerator of propagators, i.e. from p ■ 7 and not the mass 
term. The last two graphs of fig. (2) contain continuous lines of fermions with one internal 
propagator, such that the divergent part is contained in CV,- only. A is finite for these two 
graphs, and thus does not contribute to the beta function. More generally [17], to any 
order of the perturbation theory, any graph containing a open fermion line, which meets 
an even number of vertices, does not contribute to the beta functions of the model. One 
needs an even number of internal lines for the product of gamma matrices (appearing in 
p • 7) to give a diverging term with a non- vanishing trace. 

3.2 Bet a- functions 

The divergent one-loop correction to the four-fermion interactions are then given by the 
first two graphs of fig. (2), which are as follows. 

• For the flavour preserving four-fermion interaction: 

(i) Graphs with the insertion of the one- loop scalar self-energy: both flavours con- 
tribute to the fermion loop, which induces a factor -4 for the trace over Dirac indices. 
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The contribution is then, 

-4i 2 (iV2g a )% a - 4i 2 (iV2g a ) 2 (iV2g b ) 2 I bb = lQg 2 a (g 2 Jaa + gthb) ■ (36) 

(ii) Graphs with the insertion of the one-loop Yukawa interaction: only the flavour 
a plays a role, and the contribution is 

2i 2 (iV2g a ) 4 I aa = -8g% a . (37) 

The total contribution must be identified with to correction to the bare graph i(i\/2g a ) 2 , 
such that 

i5g 2 a = -4g 2 a (g 2 Jaa + 2g 2 I bb ) , (38) 
and the corresponding beta function is therefore 



9(Sg 2 a ) = _gl, 2 
5A ir 2 



A _U&_ = _*L( g * + 2 gt) . (39) 



• For the flavour-changing interaction: 

(i) Graphs with the insertion of the one- loop scalar self-energy: 

16g a g b (g 2 a I aa + g 2 I bb ) ; (40) 

(ii) Graphs with the insertion of the one- loop Yukawa interaction: 

-Agwlha - ±g a glhb ; (41) 
The total contribution must be identified with the correction to the bare graph iG 

t5G = -12g a g b (gXa + glhb) , (42) 
and the corresponding beta function is therefore 

/3 G = -3^f (g 2 a + g 2 b ). (43) 

One can conclude from this one-loop analysis that the theory is asymptotically free, since 
higer orders also diverge at most logarithmically, and cannot change the sign of the one- 
loop beta functions. Note that, when g 1 = g 2 , then /3 G = 2j3 a , as expected from the 0(2) 
symmetry. 



4 Two-loop propagator 

Since Lifshitz theories explicitly break Lorentz symmetry, space and time derivatives are 
dressed differently by quantum corrections. If one considers only one particle, or sev- 
eral particles in a given flavour multiplet, frequency and space momentum can always be 
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a,b a 



Figure 3: Two-loop contributions to the propagator. The fermion loop in the second graph 
involves a contribution from each flavour, and a factor -4 for the trace over Dirac indices. 



rescaled in such a way that the particles have the usual Lorentz-like IR dispersion relation 
(after neglecting the higher order powers of the space momentum, suppressed by M). How- 
ever if one considers several particles without flavour symmetry, then it becomes necessary 
to perform a flavour-independent rescaling of frequency and space momentum, such that 
different particles see different effective light cones. This is the case we consider here. 
The dispersion relations (4) are not modified at one loop, since the corresponding graphs 
do not depend on the external momentum. We therefore have to go to two loops to find 
the first quantum corrections, corresponding to the graphs represented on fig. (3), in terms 
of the equivalent Yukawa model (34). 

We note here that the two-loop propagator is evaluated in [18] for a scalar 4 theory, in 6 
spacial dimensions and for z — 2. This calculation is done in the massless case and in the 
absence of quadratic space derivatives. Dimensional regularization is used there, such that 
the power of the cut off does not appear explicitly in the results. The authors conclude 
that the the Lorentz-symmetry breaking terms flow to in the deep IR. 



4.1 Self energy 

The perturbative graphs on fig. (3) can be calculated with massless bare propagators, 
since the two-loop graphs contain no IR divergence. As a consequence, these graphs are 
flavour independent (besides an overall factor depending on the coupling constants), and 
they involve the integrals 

T(u u\ - . 2 f dppdp r dqpdq iN(-p)iN(-q)iN(p + q + k) 
1 °' ' ~ 1 J (2tt) 4 7 (2tt)4 D(-p)D(-q)D(p + q + k) [ ] 

dpodp f dqodq Tr[iiV (—p)iN (—q)]iN \p + q + k) 



(2tt) 4 J (2tt) 4 D(-p)D(-q)D(p + q + k) 

where the trace in J arises from the fermion loop, and the factors i 2 are for the scalar 
propagators. Taking into account the different possibilities for the self-energy of flavour a, 
we obtain 

• (iy/^ga) 4 ! for the graph without fermion loop; 
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• [(iy / 2fi'a) 4 + (iV^ga) 2 (iV^gb) 2 ]J for the graph with a fermion loop: one contribution 
for each flavour in the loop. 

We calculate these integrals in the Appendix, where we see that the only role of the fermion 
loop is to give a factor -4 from the trace over Dirac indices. We therefore have J = —47, 
and the total contribution to the momentum-dependent two-loop self energy T> a (k ,k) is 
given by 

-tZ a (k ,k) (45) 
= -452(3^ + 4^)7 

dp dp f dq dq N(—p)N(—q)N(p + q + k) 



-^g a (3g a + 4g bJ j ^ 4 j ^ D (_ p ) D (_ q ) D (p + q + k ) ■ 
The bare inverse fermion propagator is 

S^ re = k ol °-M 2 k^+--- , (46) 

— * 

where dots represent higher orders in k. We parametrize the dressed inverse propagator as 

Sussed = -ml + (1 - Y a )k ol ° - (1 - Z a )M 2 k • 7 + • • • , (47) 

such that the self energy is 

£ a (A; , k) = S^ re - S~ dr \ ssed = ml + Y a k ol ° - Z a M 2 k ■ 7 + • • • (48) 

The integrals (44) should then be expanded in the external frequency A; and momentum 
k in order to find the corrections Y ai Z a . The /^-independent mass correction m a will be 
disregarded, since the dynamical masses have already been calculated. 

4.2 Dressed dispersion relations 

From the self energy (48), the IR dispersion relation for the flavour a is 

(1 - Y a fkl =mt + M 4 (l - Z a fk 2 + • • • , (49) 

— * 

where k = \k\. If we assume that the two fermion flavours are to be coupled to other 
particles, then one needs a flavour-independent rescaling of the dispersion relation, ko — > 
M 2 ko leads then to the following product of the phase and the group velocities, v p and v g 
respectively 

v 2 a = v p v g = = 1 + 2{Y a - Z a ) + 0{k/Mf . (50) 

We calculate Y a and Z a in the Appendix, by expanding analytically the integral 7 to first 
order in ko and k ■ 7, and we find a quadratic divergence of the form (A >> M) 

Y a -Z a ^ Ak g 2 a (3g 2 a + 4g 2 )^ , k ~ -3.49 x 10~ 5 (51) 
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where a ^ b. This result shows that the present model is of limited use for the prediction 
of Lorentz-violating propagation. Indeed, with the values of A/M, gi, g 2 shown in Table 
1, the result (51) is not perturbative: one needs to absorb the quadratic divergence with 
counterterms, such that the renormalized value of Y a — Z a needs to be fixed by experimental 
data. Therefore the model cannot predict quantitative deviations from Special Relativity 
at low energies. 

If these corrections were logarithmic, one could infer from our result a cut-off-independent 
beta function for the effective maximum speed seen by the fermions, which could lead to 
"realistic" predictions on potential sub/super-luminal propagation. 

Note that the rescaling of frequency which leads to the speed squared (50) does not make 
apparent the fact that, if flavour symmetry is exactly satisfied, then the IR dispersion 
relation are relativistic. If one ignores possible interactions with other particles, one can 
further rescale 

k 2 = ff i-ri 1-^ (52) 

1 _ Zl 1 - z 2 ' K J 

which leads to the following IR dispersion relations 

~k 2 ~ m\ + (1 + 25v)k 2 
k 2 ~ rh 2 2 + (l- 25v)k 2 
A 2 

where 5v = 6k(^-^|)— . (53) 

One can see here that the Lorentz-invariant IR dispersion relations are recovered when 
g 1 — g 2 . But for g x ^ g 2 , one needs the difference \gf — g 2 \ to be proportional to M 2 / A 2 in 
order to deal with realistic phenomenology. One can take the example of the largest value 
M/A ~ 10~ 5 with g a ~ 0.58 from Table 1. The upper bound 5v < 2 x 10~ 9 given by the 
supernovae SN1987a data [19] gives then 

. 9 9| Sv M 2 /A 2 . . 

|J - 9||£ M5FTS)^ (54) 

such that flavour symmetry can be considered exact, and the corresponding fine tuning is 
not natural. 



5 Conclusion 

In the Lifshitz context, we have shown that flavour oscillations are generated dynamically, 
for massless fermions coupled via four-fermion interactions. The corresponding IR disper- 
sion relation is not relativistically invariant, as a consequence of the absence of flavour 
symmetry. An essential point is that, even if Lifshitz kinematics are hardly detectable at 
the classical level, quantum corrections lead to important effects in the IR, as flavour oscil- 
lations and modified dispersion relations. Nevertheless, we find that modifications to the 
IR dispersion relations lead to too important Lorentz-violating effects, which are obviously 
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ruled out, unless flavour symmetry is almost exactly satisfied, which corresponds to an 
unnatural fine tuning. 

We therefore suggest that a realistic Lifshitz model should have logarithmic divergences 
at most, in order to have phenomenological relevance. This is the case, for example, of 
Lifshitz- type Yukawa models [20], where one-loop corrections to the fermion dispersion 
relations are finite. Also, Lifshitz-type extensions of gauge theories, which are super- 
renormalizable in 3+1 dimensions and for z — 3, feature interesting properties [21]. An es- 
sential point is, since the gauge coupling constant has a (positive) mass dimension, fermion 
dynamical mass naturally appears in the dressed theory. If, in addition, fermion conden- 
sates break flavour symmetry, then vectors automatically become massive [22], with a 
mechanism similar to the one initially derived by Schwinger in 1+1 dimensional Quantum 
Electrodynamics, and used in Technicolor studies. Such models are to be looked at in 
future publications, in order to determine their phenomenological relevance. 

Acknowledgements We would like to thank Antonis Tsapalis for useful discussions on 
the four-fermion model, and Alessio Comisso for suggestions on numerical integrations. 



Appendix: two-loop propagator 



An expansion in the frequency k of the integrand appearing on the right-hand side of 
eq.(45) gives 



N(-p , -p)N(-q , -q)N(po + q + k ,p + q) = N(-p)N(-q)N(p + g) 
D(-po, -p)D(-q , -q)D(p + q + k 0: p + q) D(-p)D(-q)D(p + q) 
o N(-p)N(-q) 



(55) 



Ohfr, \ N(-p)N(-q)N(p + q) 2 
2fc (Po + qo)-^ rrr7 TT^T — ; — r + O(k ) , 



D(-p)D(-q)D(p + q) " uv " u ' " u/ D(-p)D(-q)D 2 (p + q) 
where (p) = (po,p), and an expansion in the spatial momentum k gives 
N(-p , -p)N(-q , -q)N(p + q ,p + q + k) N(-p)N(-q)N(p + q) 



D(-p , -p)D(-q , -q)D(p + q ,p + q + k) 
N(-p)N(-q) 



+ 



D(-p)D(-q)D(p + q) 
N(-p)N(-q)N(p + q) 
D{-p)D(-q)D*(p + q) 



D(-p)D(-q)D(p + q) 
2((p + q) ■ k)((P + Q)-l) + (M 2 + (p + q) 2 ){k ■ 7) 



(56) 



2((p + q)-k) 



M 



(p + q) 2 M 2 + 3(p + q) 2 + O(k 2 ) . 



The first term in the /^-expansion leads to the integral 

N(-p)N(-q) 



jp 



(57) 

? iy\— V) iy\—K l) iyyy t qj 

/ i rO — 0\ / 1 rO — 0\ — » — » 

= (hi ) I 

Jv,< 



P , q D(-p)D(-q)D(p + q) 

p q -(M 2 +f)(M 2 + f)p-q 
l P , q D(-p)D(-q)D(p + q) 
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and, because of the symmetry p •<->■ q, the second integral vanishes. The following rescaling 
Po = M 3 u , q = M 3 v , p = Mu , q = Mv, (58) 
together with a Wick rotation on Wo^o finally leads to 

(U ^ f N(-p)N(-q) 

{k °' r) LD(- P )D(- q )D( P + q ) (59) 

on f d^UE d 4 V E M 4 f 4 + (1 + U 2 )(l + V*)u • V . , 

lfc ° 7 j i (27r) 4 (27r) 4 L> E (^)L> £ (^)L> £ ( WB + ^) ' 1 j 

where D E {u E ) — u\ + (1 + u 2 ) 2 u 2 . 

The second term in the A; -expansion gives 

, /■ , , x N(-p)N(-q)N(p + q) 

/■ go/ g • 7 (p + g ■ 7 (M 2 + g 2 )[M 2 + (p + g) 2 ] - (p*+g) 

-2A; y^Cpo + go) D(-p)D(- g )D'(p + g ) ' 

where, by symmetry, the terms proportional to 7 lead to a vanishing integral. After the 
rescaling (58) and a Wick rotation, we then obtain 

, f ( , , N(-p)N(-q)N(p + q) 

o(b °\ f ^ UE d4yE ( -i- \ 2 um + + ^ 2 ^ 1 + : 

" Afc ° 7 j 7 (2vr) 4 (2tt) 4 ^ + ^ + v E ) ' 

The term proportional to /c 7° is then 

/, On f <^UE_(^Ve_ M 4 ^4 + (1+W 2 )(1+^)M-^ /_ 2(tt 4 + ?7 4 ) 2 \ 

l/C ° 7 j 7 (2tt) 4 (2tt) 4 L> £ (^)L> E K)L> E (^ + ^ + D E (u E + v E )J ' ^ 

— * 

For the first term in the /c-expansion, we use the identity 

/ f(p,q)p-k (p + q)-l = ^- [ f(p,q) p-(p + q) , (64) 

Jp,q 6 Jp,q 

where f(p, q) depends on (p) 2 , (q) 2 and p • q only. The rescaling (58) and a Wick rotation 
then lead to the integral 

- l o ( -Srs' + ,) i 2<(?+ 91 ■ * )(<?+ *> ■ + (m2 + <?+ « 2)( * ■ ^ 

f d A u E d 4 v E / 5_ , 2 \ -u 4 w 4 + (1 + w 2 )(l + v 2 )^ ■ v . . 

= M 2 (k-^) - — =-- — =j [l + -(u + v) 2 -— — r— — rr— -. (65) 

1 1J J (27r) 4 (27r) 4 V 3 l >) D E (u E )D E (v E )D E (u E + v E ) 1 ; 
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The second term in the ^-expansion leads to the integral 
N(-p)N(-q)N(p + q) 



i 



P , q D(-p)D{-q)D*(p + q) 



((p + q)-k) M 2 + (p + q} 2 M + 3(p + q) 



= -2M 2 



u v — (1 + u 2 )(l + v 2 )^ ■ V 



(1 + (u + vf)(u + v) .7 



L,v D e (ue)D e (ve)D 2 e (u e + v E ) 

x(u + v)-k [1 + (u + v) 2 ][l + 3{u + v) 2 } , (66) 

where, by symmetry, the term not proportional to 7 vanishes. Using the identity (64), a 
Wick rotation then leads to 

N(-p)N(-q)N(p + q) _ ^ , 2 + , , 2 + 3 



(67) 



Pl ,I>(-p)£>(-g)£> 2 (p + g) 
2 M 2 (k ^\ I ^ 4 u E d 4 v E u 4 v 4 + (I + u 2 )(l + v 2 )^ ■ v 
"3 ( ' 7j 7 (2tt) 4 (2tt) 4 D e (ue)D e (ve)D e (ue + v E ) 
x (« + t;) 2 [l + (« + {;) 2 ] 2 [l + 3(« + ^) 2 ] . 

The term proportional to (k ■ 7) is then 

d 4 u E d 4 v E U4V4 + (1 + m 2 )(1 + v 2 )u ■ v 



M 2 (k-i) 



(2tt) 4 (2tt) 4 D e (u e )D e (v e )D e (u e + v E ) 

x ( 1 + ^+v) 2 --(u + v) 2[ -^±^- 

V 3 l ; 3 l ; D E (u E + v E ) 



\212 



[1 + 3(m + {7) 2 



(68) 



Finally, from eqs. (63,68), the quantum corrections to the IR dispersion relation are deter- 
mined by 

Z a = 4g 2 a (3g 2 a + 4g 2 ) 



Yn 



d 4 u E d 4 v E 



(2tt) 4 (2tt) / 



Int , 



(69) 



where the integrand is 

1 U4V4 + (1 + u 2 )(l + v 2 )u ■ V 



Int 



3 D E (u E )D E (v E )D 2 E (u E + v E ) 



(70) 



x 



(« 4 + ^) 2 (6 + 5(u + t;) 2 ) - (iT + t/) 2 (l + (u + v) 2 f{2 + (u + v) 2 ) 



Note that in the Lorentz-symmetric case, higher orders in u, v are absent and Y a = Z a . 
The integral (69) is evaluated as follows. 

We can first perform the exact integration over -u 4 , t> 4 , using the Feynman parametrization. 
This introduces two new variables of integration, but which lie in a compact domain of 
integration: 

1 



D E {u E )D E (v E )D 2 E {u E + v E ) 

r-l rl—x 

— 6 / dx I dy 



1 



x-y 



[xD E (u E ) + yD E (v E ) + (1 - x - y)D E (u E + v E )\ 



4 ' 
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We then introduce the variables a, b, such that 



■u 4 = s(a + b) and t> 4 = t(a — b) , with s = y/\ — x and t — yl — y , 
to obtain 

d«4 f 1 , f 1 ^ x , 12 da d& sta 

D e (u e )D e (ve)D 2 e (u e + u B ) 

where 



rl rl—x 

dx dy 
Jo Jo 



[2st(st + a)a 2 + 2st(st - a)b 2 + D} 4 



D = x(l + u 2 ) V + y(l + v 2 ) V + (7(1 + E) 2 E 
£ = (-u + {7) 2 , a" = 1 — rr — y . 



We then write, with O<p<oo,O<0<27r 



to obtain 



2st(st + cr) a = pcos0 and J2st(st — a) b = psin0 



y d-U4 y di>4 Int 

2 io ^/o dy 7W^Jo pdp L d(t) WTW 



X 



X 



p 2 / cos 2 sin 2 (f) 



+ (1 + M 2 )(1 +W 2 )m-{T 



2 y s£ + cr st — a 
p 2 /(s + t) 2 cos 2 , (s-t) 2 sin 2 
2st V 



st + cr 



+ 



/■l rl—x 

= 2ix \ dx \ dy 
Jo Jo 

= — [ dx [ dy- 
6 Jo Jo 



a 



y/s 2 t 2 - a 2 Jo 



a 



st — a 



°° , Ap A + Bp 2 + C 

pdp [p 2 + DY 
2A B 2C 



(6 + 5S) -S(l + S) 2 (2 + E) 



where 



Vs 2 t 2 - a 2 V D D 2 D 3 
6 + 5S 2s 2 t 2 + Aa 2 - 3a(s 2 + t 2 



( S 2 t 2 _ a 2)2 



s 2 + t l - 2a 



s 2 t 2 - a 2 



B = (1 + u 2 )(l + v 2 )u ■ v(Q + 5E) 
C = -2(1+m 2 )(1 + w 2 )m-{;S(1 + S) 2 (2 + S) . 

We then define u ■ v = uv cos and 



+ S(l + S) 2 (2 + S) 



a 



s 2 t 2 - a 2 



u = rcosa , f = rsino; , with < r < oo and < a < it/2 , 
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and the final integral is 

= 1 Z 1 dx t~ X dy [ A/M dr r d9 da ^a±^ (*A + A + ^ 

which is quadratically divergent, as F(z) ~ /tz 2 when 2—7-00. We then find via numerical 
integration 

f 1 dF) 

k = lim s — I 80 

\ 2z dz I v 7 

f r 4 /-i r 1 -^ , ,„ r /2 r * sin 2 (2a)sinfl (2A B 2C\\ 
= hm { — - — - dx dy d0 da =j= = — + — + 



r— >oo 



3x2% 5 i Jo w jo jo Vs^ 2 - 2 \ D D 2 J J 

3.49 x KT 5 , (to a 1% accuracy) . (81) 
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